1 



DUKE-TH-02-219 



(N 
O 
O 



O 

(N . 

(N ■ 
> 

O . 
(N ■ 

O ■ 
cn : 
o . 

o : 

Oh; 

^ ■ 



X 



Unexpected Results in the Chiral Limit with Staggered Fermions 

Shailesh Chandrasekharan ^ 

''Department of Physics, Duke University, Durham, NC 27708-0305, USA 

A cluster algorithm is constructed and applied to study the chiral limit of the strongly coupled lattice Schwinger model involving 
staggered fermions. The algorithm is based on a novel loop representation of the model. Finite size scaling of the chiral susceptibility 
based on data from lattices of size up to 64 x 64 indicates the absence of long range correlations at strong couplings. Assuming 
that there is no phase transition at a weaker coupling, the results imply that all mesons acquire a mass at non-zero lattice spacings. 
Although this does not violate any known physics, it is surprising since typically one expects a single pion to remain massless at 
non-zero lattice spacings in the staggered fermion formulation. 



1. Motivation 

Consider a two dimensional U{1) lattice gauge theory 
with staggered fermions. This model has been studied 
extensively on the lattice for over two decades In 
the limit when the lattice gauge coupling goes to zero, 
the model is believed to describe the continuum two fla- 
vor Schwinger model; the fermions are confined and the 
low energy physical particles are mesons. Further, in the 
chiral limit three massless and one massive pseudo-scalar 
mesons emerge like in QCD §]. However, since in two di- 
mensions continuous chiral symmetries remain unbroken, 
the chiral condensate vanishes in the chiral limit. The 
actual prediction is 



,1/3 



(1) 



as discussed in |^ . This means that the chiral susceptibil- 
ity X = ^ ('0V') diverges in the chiral limit reflecting the 
presence of massless excitations. This is consistent with 



the fact that 



m 



2/3 



, i.e., the pion mass vanishes at 



TO = 0. If one computes the susceptibility at to = in a 
finite box of size L x L, one expects 

X - L, (2) 

i.e., it diverges linearly with L pO[ |. 

When the lattice gauge coupling is non-zero, lattice ar- 
tifacts in the staggered fermion formulation break the chi- 
ral symmetry of the two flavor Schwinger model explicitly 
down to a U{1) subgroup. As far as we know, no one 
has completely analyzed the effects of lattice artifacts on 
the particle spectrum. Do massless particles remain in 
the spectrum? The conventional wisdom from higher di- 
mensions is that there should be one massless pion, since 
the remnant t7(l) chiral symmetry is expected to break 
spontaneously at strong couplings. Of course the Mermin- 



Wagner theorem forbids the breaking of a continuous sym- 
metry in two dimensions, although non-interacting mass- 
less bosons can emerge. On the other hand a U{1) (or 
equivalently an 0(2)) symmetry is special. In such a case 
bosons can interact through topological excitations as dis- 
covered by Kosterlitz and Thouless ||l^. Thus, the lattice 
model could either be in a massive or a massless phase 
depending on the effective couplings of the low energy ef- 
fective theory describing the bosonic excitations. In the 
massless phase there are predictions for the behavior of 
the chiral condensate and the susceptibility based on uni- 
versality p^. One expects 



J7/(4-?)) 



(3) 



where < ry < 0.25. Again the susceptibility x diverges 
in the chiral limit reflecting the presence of massless exci- 
tations. In this case the finite size scaling formula for the 
susceptibility is given by 



X 



(4) 



It is amusing that if we set r/ = 1 in eqs. and we 
recover eqs. (|l|) and (||). 

Most results from earlier work on the lattice Schwinger 
model favor the view point that there is one massless pion 
at finite lattice spacings. This is based on the observation 
that the pion mass decreases with the fermion mass like 
in the continuum Schwinger model. However, on closer 
examination the lattice model shows deviations, which 
become larger at stronger couplings as one might expect 
On the other hand, no one has ever found scahng 
relations expected from universality in a massless phase 
of an 0(2) model in two dimensions. In particular no one 
has been able to confirm relations similar to eqs. (^) or 
(U). Thus, inspite of the large amount of literature on the 
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subject, questions related to the chiral limit of the lat- 
tice Schwinger model with staggered fermions still remain 
unanswered. 

The essential difficulty is the absence of a reliable nu- 
merical approach to study interacting fermionic field the- 
ories close to the chiral limit. For this reason most pre- 
vious studies alluded to above, have focused on calcula- 
tions away from the chiral limit and have used extrapo- 
lation techniques to predict the chiral limit. As we will 
see this can be misleading. In the last few years fermion 
algorithms based on cluster updates have emerged, which 
do not suffer from problems that conventional algorithms 
face near the chiral limit |13|. One can now work directly 
in the chiral limit using the new approach, a luxury not 
available with earlier methods. Recently, a fermion clus- 
ter algorithm has confirmed the scaling predictions near 
a Kosterlitz-Thouless transition in a Hubbard type model 
with unmatched precision Interestingly, these new 

ideas can also be applied to study the chiral limit of the 
lattice Schwinger model at strong couplings The al- 
gorithm is based on a novel loop representation of the 
model. Although, it may be possible to extend the method 
to weaker couplings, in this article we concentrate on the 
strong coupling limit. We focus on the question of whether 
there are massless excitations at strong couplings by look- 
ing for a divergence in the chiral susceptibility as predicted 
by eq.(i. 

2. The Model 

The two dimensional U{1) lattice gauge theory with 
staggered fermions is described by the action 



[] (1 -f V.V'x^.H-^^: 



(6) 



which can be written as a sum over weights of configu- 
rations of monomers and dimers [p^ . Each configuration 
consists of Tlx — 0,1 monomers on site x and bx,f_i = 0,1 
dimers on the bond connecting x and x + /I. In order for 
the Grassmann integration to give a non-zero result we 
need 



nx + bx^i + hx,2 + b^_i ^ + 6^_2 



1. 



(7) 



x-p,) 



Figure 1. A typical monomer- dimer configuration satis- 
fying eq.^. The weight of the configuration is m^. 



'"Y^x^x, 



(5) 



where x = {xi,X2) represents a lattice site on an L x i 
lattice, fijiy = 1,2 represent the two directions and fi, v are 
the unit vector in the /i and the v direction respectively. 
The phase factors r}i{x) = 1 and 772(2;) — (—1)^^ are the 
staggered fermion phase factors. The gauge fields Ufi{x) 
are elements of U{1) phase factors and the fermion fields 
i^ix) and "0(2^) are Grassmann numbers. 

In the strong coupling limit {g — 00) it is possible to 
integrate over the gauge fields and obtain 



z(m) = / [dVlMV'] n 



1 



Assuming this constraint the partition function can be 

written as 



Z{m) 



[nM 



exp 



log(m) 



(8) 



A typical configuration is shown in figure |l|. 

When m = no monomers are allowed and the partition 
function is given by the number of closely-packed-dimer 
(CPD) configurations on the lattice. Such configurations 
are interesting even in condensed matter physics and have 
played an important role in the study of the 2-d Ising 
model ||l^. The chiral symmetry of staggered fermions 
is manifest at finite volumes through the fact that the 
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chiral condensate vanishes since it is impossible to find 
a CPD configuration with just one monomer. The chiral 
susceptibility on the other hand is non-zero and is a useful 
observable in the chiral limit. It is equal to the total num- 
ber of CPD configurations with two monomers, where one 
of the monomers is constrained to be at a fixed position, 
divided by the partition function. 



rules of one such extension. If we ignore monomers each 
shaded plaquette of the configuration of Fig. |l| carries one 
of the seven plaquette configurations given on the left side 
of equations in Fig. ||. On the other hand, the right side of 
the equations represent the new configurations and their 
weights. Figure |^ gives an example of a loop configura- 
tion. In the absence of monomers it is easy to check that 



3. Loop Representation 

It is easy to construct a local Metropolis algorithm for 
the monomer-dimer model when m ^ 0. The algorithm is 
based on an update which either breaks a dimer into two 
monomers or vice- versa. This algorithm works reasonably 
well for TO > 0.01. However, the algorithm fails in the 
chiral limit since no monomers are allowed when m = 0. 
In fact in the chiral limit it is easy to find configurations 
where simple local updates do not lead to another allowed 
configuration. Perhaps for this reason, as far as we know, 
no one has successfully studied the chiral limit. 



Figure 2. Rules for extending the CPD configurations to 
include additional bond variables. 

To construct an algorithm that is applicable in the chi- 
ral limit we extend CPD configurations to configurations 
of loops made up of bonds which include the original or 
"filled" dimers (represented here by "solid" bonds) and 
"empty" dimers (represented by "dashed" bonds) such 
that the partition function can be expressed as a sum over 
weights of new loop configurations. Figure shows the 




Figure 3. An example of a loop configuration with two 
monomers. 

all constraints are satisfied if each loop is made up of a re- 
peating sequence of filled and empty dimers. This means 
that there are two allowed configurations associated with 
each loop. The usefulness of the loop variable is that a 
dimer system can be updated by "flipping" a loop where 
filled dimers are emptied and vice versa. The Metropolis 
acceptance of such a loop flip is found to be reasonable al- 
though large loops are not flipped as often as small loops. 
When monomers are allowed each loop consists of an even 
number of them. Further, when monomers are present 
in a loop then it contains a unique pattern of filled and 
empty dimers and a flip is not allowed. However, close to 
the chiral limit loops containing monomers are negligible. 

The algorithm to measure the susceptibility in the chi- 
ral limit is straight forward in the loop representation. 
Typically, we choose a point at random on the lattice and 
traverse the loop it is on, starting in the direction of a filled 
bond. As one goes around the loop, the bonds are flipped 
and the change in the weight of the configuration is noted. 
Interestingly, every time a filled bond is emptied one gets a 
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configuration that contributes to the susceptibility. Since 
its weight relative to the original configuration is known 
at the time of the flip, it is recorded as a part of the mea- 
surement. When the whole loop is flipped one knows the 
change in the weight of the configuration and a Metropolis 
accept-reject step can be performed. If the configuration 
is accepted then one gets a new configuration. Otherwise 
the loop is flipped back to the original configuration. In 
any case, the sum over all the weight changes recorded 
while flipping the filled bonds in the loop divided by one 
plus the total weight change due to the loop fiip is taken 
as a contribution to the susceptibility during that update. 
It is possible to show that this algorithm is crgodic. 

4. Results 

Since the remnant U{1) chiral symmetry of staggered 
fermions can break spontaneously in higher dimensions, 
the chiral condensate can approach a constant in the chiral 
limit. However, in two dimensions a continuous symmetry 
cannot break and it is almost guaranteed that 



(9) 



with < 5 < 1. Comparing with cqs. (|l|) and (^) we see 
that S = 1/3 in the continuum two fiavor Schwinger model 
and d < 1/15 is expected from 0{2) universality in two 
dimensions. It is important to remember that the U{1) 
chiral symmetry of staggered fermions is non-anomalous. 
Thus, a non-zero expectation value of (tpip) in the chiral 
limit would definitely imply a spontaneous breaking of the 
U{1) symmetry. 
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Figure 4. Chiral Condensate as a function of inverse 
volume at various masses. 



Using the local Metropolis algorithm we calculated the 
chiral condensate for five masses in the region 0.01 < m < 
0.1 for L = 16,24,32 and 64. Below m = 0.01 the algo- 
rithm slows down and is not very efficient. Our results 
are shown in figure ^. The solid lines are linear fits of the 
data to the finite size scaling formula (ipip) = A + B/L'^ 
at a fixed value of m. The constant A then yields the 
condensate at infinite volume for a fixed mass. Figure ^ 
shows the condensate as a function of the mass for L = 16 
and for infinite L obtained from fits shown in figure ^. At 
L = 16 we see that the condensate goes linearly to zero as 
expected due to finite volume effects]]. On the other hand 
the infinite volume results fit beautifully to a power law 
of the form (ipip) = 0.455(6)m°-^°^^^). Since we are in the 
strong coupling limit it is not surprising that the power 
does not match the continuum two flavor result of eq. (|^) . 
However, it neither matches the predictions of the two di- 
mensional 0(2) model based on universality (see eq.P)). 
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Figure 5. Chiral condensate as a function of mass. 

Before attempting to understand the unexpected value 
of the power it is useful to confirm that the power law 
behavior is valid all the way to the chiral limit. A power 
law for the chiral condensate implies that the chiral sus- 
ceptibility will diverge at m = in the thermodynamic 
limit. This divergence typically manifests itself in a finite 
size scaling of the form 



X 



_d_ 

dm 



L 



1 



(10) 



^ In a finite box the partition function Z(m) is a polynomial of the 
variable 'm? . The condensate is proportional to the first derivative 
of Z(m) with respect to m and hence vanishes linearly with m. 
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consistent with eqs. ^ and (||). This behavior is diffi- 
cult to study with conventional algorithms. However, the 
loop cluster algorithm discussed in the earlier section is 
designed exactly for this purpose. In figure ^ we plot our 
results for the chiral susceptibility obtained with the new 
algorithm at m = as a function of the lattice size. Notice 




Lattice Size (L) 



Figure 6. Finite size scaling of the chiral susceptibility. 

that for L = 16 the value of the susceptibility is consistent 
with 15, the slope of the condensate at m = shown in 
figure ||. 

As a function of L the susceptibility behaves differently 
in different regions. The data is described well by a power 
law in the region 8 < L < 20 with 7 = 1.56, while a 
straight line fits the data in the region 20 < L < 40. For 
L > 40 the susceptibility appears to be approaching a 
constant. If we compare the power law behavior of the 
susceptibility in the region 8 < L < 20 with eq. we 
find that 77 ~ 0.44 which is inconsistent with the pre- 
dictions of 0(2) universality in which one would expect 
< r? < 0.25. However, it explains the strange power 
law behavior of the infinite volume chiral condensate we 
encountered in figure |^. If we substitute 77 = 0.44 in eq. 
(H) we find that (ip^) mP'^^, which is quite close to the 
L = 00 results of figure || in the region 0.01 < m < 0.1. 
Since the power law behavior of the susceptibility does not 
extend to larger volumes, it is clear that the power law be- 
havior of the infinite volume chiral condensate shown in 
figure H cannot hold all the way to the chiral limit. In fact 
the susceptibility at m = saturates at large volumes. 
Equivalently, the first derivative of the chiral condensate 
with respect to the mass at m = reaches a constant as 



the volume becomes largcQ This means the infinite vol- 
ume chiral condensate will also vanish linearly at m = 0. 
Typically, divergence of a susceptibility such as the one 
defined in eq. (|l^) , is related to the presence of massless 
particles in the theory. Our results indicate that there are 
no such particles. 

5. Discussion 

The results that we have obtained are somewhat sur- 
prising. The conventional wisdom is that strong couplings 
break chiral symmetry at least in three or more dimen- 
sions. This suggests that the low energy effective model 
that describes the U{1) chiral dynamics of the lattice 
Schwinger model with staggered fermions is most likely in 
the low temperature (massless) phase at strong couplings. 
This reasoning appears to be in contradiction with our 
findings since we do not find the expected divergence de- 
scribed by eq. which in turn implies that there are no 
massless particles at strong couplings. 

Is it possible we have missed something? Of course, we 
rely heavily on our algorithm. It is rather new and may 
have weaknesses like long auto-correlation times that we 
have not yet appreciated. In order to alleviate such fears 
we have compared our Monte Carlo results with exact cal- 
culations of the susceptibility at m = on small lattices. 
Table |^ shows this comparison. Clearly, we can reproduce 

Table 1 

Chiral Susceptibility: algorithm vs. exact results. 



Lattice Size 


Exact 


Algorithm 


4x4 


1.70588235... 


1.7059(1) 


6x6 


3.33640880... 


3.3364(1) 


6x8 


4.07961565... 


4.0796(2) 


8x8 


5.27221660... 


5.2722(2) 



the exact results with great precision at least on small lat- 
tices. We have not seen any pathologies in the simulations 
at larger volumes except for the fact that the fluctuations 
increase which require us to increase the statistics propor- 
tionally. Based on this we are prejudiced to believe in our 
estimate of the errors within a factor of two or three. 

Assuming our results are correct, we find that at strong 
couplings the lattice model is in the high temperature 
(massive) phase of the low energy effective model. We 
know from experience in higher dimensions that chiral 
symmetry breaking effects become weaker at smaller cou- 
plings. This, coupled with the fact that there is no decon- 
finement transition in two dimensions, suggests that it is 

^Although we cannot rule out a mild divergence of the susceptibility 
in figure H, we think it is unlikely. 
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unlikely that there is a phase transition to a massless phase 
at weaker couplings. Our results then imply that there are 
no massless excitations for all g 7^ (or equivalently at fi- 
nite lattice spacings); the non-singlet pseudo-scalar meson 
(which would be expected to become the Goldstone bo- 
son in higher dimensions ) is actually massive at all non- 
zero couplings. This conclusion may come as a surprise 
to some, but it does not contradict any known physics 
as far as we can tell. In particular it does not contra- 
dict our expectations that in the zero coupling limit one 
must recover the continuum two flavor Schwinger model. 
It is well known that two other non-singlet pseudo-scalar 
mesons acquire a mass due to lattice artifacts in the stag- 
gered fermion formulation. Our data is suggesting that all 
mesons become massive at non-zero couplings contrary to 
expectations from higher dimensions where the Goldstone 
boson will remain massless at finite lattice spacings. In 
the lattice Schwinger model this boson is still perhaps the 
lightest. The gauge coupling controls the lattice spacing 
in the model; typically one uses the mass of the iso-singlet 
meson to set the lattice spacing. This means that in order 
to recover the correct continuum limit the singlet and the 
non-singlet meson masses must become smaller in lattice 
units as one approaches weaker couplings but their ratio 
must diverge. 

The current work can be extended in many directions. 
There are results in the large Nc and large d limits at 
strong couplings which show that chiral symmetry is in- 
deed broken [|l8|;^. Here we have shown that at A'c = 1 
and d — 2 chiral symmetry remains unbroken and the the- 
ory is in the massive phase. It would be interesting to find 
the value of Nc for d = 2 at which the theory will move 
into a phase with massless excitations. There is a lot of 
evidence from lattice simulations and chiral extrapolations 
that shows lattice QCD with staggered fermions to be in 
the chirally broken phase. However, as we have seen in 
the present context, chiral extrapolations can sometimes 
be misleading. Fortunately, cluster algorithms of the type 
used here can be developed for any Nc and d so that these 
questions can be answered directly in the chiral limit. One 
can start with the results presented in ||2^]. It would be 
useful to confirm that chiral symmetry is indeed broken at 
Nc — 3 and d — A. On the algorithmic side, it is exciting 
that we have a new method to explore the chiral limit in 
certain gauge theories. Loop variables play an important 
role in this method. Such variables have already been dis- 
covered in a limited class of models |^,^ . Perhaps it is 
possible to extend them to other interesting models. 
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